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ABSTRACT 


Allocation of resources has become a classical 
problem in optimization by mathematical programming. 
In the field of military applications attack aircraft 
assignment has been treated widely by deterministic 
and/or linear models. However, destruction of a target 
is no certainty nor is damage inflicted on a target 
linear with respect to the number of weapons delivered 
on it. Recent extensions in the field of nonlinear 
programming in conjunction with the widespread use of 
electronic digital computers permit a more realistic 
approach to this problem. This paper formulates a 
stochastic nonlinear model for assigning a force of 
attack aircraft on a single sortie against fixed location 
targets. The number of aircraft alive at weapon release 
on any pass of a series against a given target is treated 
as a random variable. The total value of damage to all 
targets is taken as the measure of effectiveness and a 
particular form of the objective function derived. The 
parameters of the model and the form of the constraint 


equations are also discussed. 


Li 





TABLE OF CONTENTS 

CHAPTER PAGE 
Tie PROB ..... . 5 cis o vis eiste 2 ieee ee Boe aL 

WOK ET OUMG,..o oo Bio. § 0s 0u Shore OME Miss buouobeuicns GE RE clcll Bicusls 

ASSUMP LOMS. .c 6.000 « 6:6 eee Gere one cme 10e o oe eens orn 

cee” OW etter sess a hia stoic « ¢ el Sav0c) o stebake aieiaioicie is a) tiels 

II. THE OBJECTIVE FUNCTION AND CONSTRAINT EQUATIONS . 


The Objective Function ©oooeoedcoeoe0000000006600 00 


oO wn wn WY NN FH 


Tne ConstraintvEQuaweions 2 crac csose clue are mnenera eieuenwiers 

Nonlinear Programs> .q.ckW Biante} h. >t ee eee Lt 

Til, "THE PROBABILITY DISTRIBUTIONS .5o.0-sccccss ue 12 

ASstmptlons . o cesses tu wtetents ob boo cso co ce 6 cee, Le 

Ween COWS . sMierepe elles 0 co c.0 0 oo pions cee cere 
Distribution of the Number of Aircraft Alive 

at Release on the First Pass ....:.scces mame 19 
Distribution of the Number of Aircraft Alive 

at Release on Subsequent PasseS ..cccccccecccos LY 

Distribution of the Total Number of Passes 

Against a Target “ssldeciintvitalslaisle come is cc o cl cmuiteme 

IV. “ERE PARAMETERS .... «cictelaieieietcteli pic oo 0c «co clei stem am 

MOY Get= VaMUe << ¥ ac c-cishete cvs o TNTe puenele Sus ole Se ee Re 

Fiiel and Ordnance Leadings 4..4 ctccececccecote 26 

Geher*Paraneters Myiewlc eles so Cals cc clctee eye ew ERED 

V. CONCLUSIONS AND RECOMMENDATIONS ...ccccccccccccce 29 


Conclusions @¢6000C0CO00 80 GCOC8KHHHCFHFHFFFFA826 FGF 0000 06 6 29 


111 





CHAPTER PAGE 
Recommendations ®eee4eeseees+e2eeoeeeee48ekreeseeneneeepeneeeoeeee &€ & @ @ 30 
BIBLIOGRAPHY @eeeeeseee21eee7een1eosvseeseeeeseee*#en,esdseeeoeeoeee7ene+eeoeeee?eee#eee@8e@ 86 @ OL 
APPENDIX A. Derivation of the Probability Distribution 
4 
of W, (N, ) @e@eeeeeoeseeeeoeoeeseoees7eseeeoe7ee1e#e#80 28 @ 3 


APPENDIX B. Program ORDLOAD @®e@eeeoeee7eeecesedseeeeee4eeseee#eeeeseeee 40 


iv 





TABLE 


If. 


Lor. 
IV. 


LIST OF TABLES 


PAGE 


Probability Distribution of X 16 


191 Ng g)eccrccores 


Probability Distribution of X (Ny 5) Given 


1 jk 

My 503 (Ny 5) = BE lo) sc oo tnteve’s oe 5 ol Le. 7 
Probability Distribution of X,,,(N,,), k #1 .. 18 
Probability Distribution of W, (N, ) eo ccccseccee 19 





LIST OF SYMBOLS 


SYMBOL DEFINITION 


D, @eeoeoeoetet @eeegsgeneo1eneeeeeeee 8 8 @ @eeoeneeteee ee & @ 
faco @eeeeeeeeeeeeneeeeeeeeeoeneeeeeeeeeee @ 
Gift ) @eeeeoeeeeeseeeoeeeeeeeseeenesesoneeneeneeee @ 
Ny ®eoeoeeeoeeeeeeeeeeweeeeeeweeeneeee7eelmtmhlUmcetmlUlcOlUlcOlmlUcPFlhlUcFlUhF8 


Ny @eeees @eeeeoeeoeeoeeeeeneneeoeeeneneneeeecmhm6UmhmUMCOMUCOOUhFMDUMO 


Pay @eeee3@#eeoeeen@eeeeeeeeenweoemlmetmlcectemlcltlUcrtrlcecCUlcrcOlcetecUlcelcemUcOmUlchOmUclhUrFhU} 
PA, y(n) eoeoeeeoeveeveoeveeeeeoev eee eee eee 8 Oe 
PQ, ,(n) eeveeneeveeveeoevoeoe eeueeonoeeneneee ee 0 e de @ 
PE, ,(N,n) e@eeeoeveveevevoeneoeveoovoeoeoeeoeoenonoenee eee @ 
PK Ce oe ee ee ee ee ee ee oe ee ee ee ee ee oe ee ee ee ee de oe ee ee et et er ee ee ee oe e@¢°¢@ 
PR, 4, (1) C@oeeoeveveeeevneeneneeneneeeeneee ee ee @ 
PIy Ce oe ee ee oe oe ee oe ee ee ee Oe oe ee oe ee ee ee ee oe oe ee ee ee ee ee er oe oe 
R @#eoeoe7e+ee2eeo@0e0eeses0e8ece0ee@e0@e@eemUmemUmUmlUCcCOUcCOCcCoee@elemUmetmlceclelUCccetclUctroeclceétclUctlh8 
Ry @oeoeoceeed@ee@eg0e@8hethlUcOlUCcCOlhUC FlC€P ®ee88 808 @ @ @®eoeoees8cecseesvse @®@ 
r @eoeoeeoese#eeee7eeseee4e01eeee#dc&edceeodse #eeo#eee¢6 0 
Te) Geaesermers 5 5 oowrarens ce errs Bas oor € 
¥, bes @#eeeseeeesesve7e2#@e#eee#esee?eee?26e?se @®9e0ese37e#e88 @ ee Clo 
Wy» W, (N, ) eoeeoeeeoeeeeeeeene eee eee ee oe 
Xs yn? Xa jac (Ny 4) Pe ore ee 


Y, @eoeoeeoeeoeee@ ®eeee @ ®©@ee8e0 0 9@e@ee 0 0088 @ 0 @ 


vi 


sGNGES...... 
D6.4 @#eee%esese?e?ees se 


$6.5 5s 7 


D2.10 @oee#eee8e86 @ 
D6.1 @e@eeeee8e080 86 


D2.3 @oeeeeee 6 @ 
D262 @#eeo7eeee#28s8e 


D6 SWORL2” . ccc 
Mo ere 
Dil ov. ee 
oe ee. 
BONS, © een 
BOs nee. 
Le 
re ee 
oe ae 
sas. oe 
ro 
Eee ccoseesen® 
D2.7 Geerccces 


D2.8 @oee3eeee#d @ 


35 
38 


38 
15, 
12 
2 
a2 


a 


BY. 


13 
7 





CHAPTER I 
THE PROBLEM 
I, BACKGROUND 


Allocation of resources has become a classical 
problem in optimization by mathematical programming. 
Linear programming has been used extensively to provide 
a solution. It is valid when the relationship between 
variables is linear, however, linearity is not always 
the case. For example, an individual may find the use- 
fulness derived from owning five automobiles not nec- 
essarily five times that derived from owning just one. 
In the past, for cases such as this, it was customary 
for one to either be satisfied with a linear approxi- 
mation or to search by other means for a Solution. The 
extension of nonlinear programming over the past five 
years to its present state provides the analyst with a 
mathematical programming technique for treating the 
allocation problem with certain types of nonlinearity 
in a direct rather than approximate fashion. Dorn fs], 
Graves and Wolfe [17] » Hadley [25] » and Wolfe [13] have 
published comprehensive works on the status of nonlinear 


programming. 





A particular allocation problem in the field of 
military applications that has received a great deal of 
attention is that of attack aircraft assignment. Numerous 
studies and papers have been published on the subject, 
and several which have been found to be of particular 
interest are? Limited War Qperations(U) [10], Tactteal 
Air Warfare Force Allocation(U) [2a], Dresher [4]. and 
Haering [29] . 


II. ASSUMPTIONS 


It was through participation in the Center for Naval 
Analyses Study. Tactical Air Warfare, 1964, that the 
necessity for a formulation such as developed in this 
paper became apparent. The allocation of attack aircraft 
to fixed location targets is treated under the following 
broad assumptions? 

(Al.1)+ An optimal allocation is one in which the 
total value of damage to all targets is a maximum, subject 
to the constraints imposed by the availability of aircraft, 
fuel, and weapons. 

(Al.2) The damage inflicted on a target is nonlinear 
with respect to the number of aircraft passes made 
against it. 


me _-_.) identifies assumptions, 
(D _._) identifies definitions. 





(Al1.3) There is a probabilistic interpretation 
of the number of aircraft alive on a particular pass 
against a given target. That is, given ten aircraft 
assigned to attack a target, there is a probability 
that all ten aircraft will not be alive on any parti- 


cular pass against that target. 
III. THE FORMULATION 


Under these assumptions the allocation of attack 
aircraft was formulated as a stochastic nonlinear 
mathematical program with its attendant objective function 
and constraint equations. The objective function in this 
case is a nonlinear function of the number of attack air- 
craft alive on each pass against each target. Its 
particular form is derived and a discussion of the con- 
straint equations given in Chapter II. 

The number of passes against a target is treated as 
a random variable in the following ways the number of 
aircraft alive at release on a particular pass against 
a target is a random variable, therefore the number of 
passes against a target is a random variable which is 
the sum of the number of aircraft alive at release on 
each pass. Their probability distributions are derived 
in Chapter III and Appendix A. 

Chapter IV is devoted to a discussion of the parameters 


of the model. These include target value, fuel required, 





racks available for ordnance, the probabilities of 
detection, engagement, kill, acquisition, and pass 
survival. 

Conclusions and recommendations comprise the final 


chapter of this paper. 








CHAPTER II 
THE OBJECTIVE FUNCTION AND CONSTRAINT EQUATIONS 
I, THE OBJECTIVE FUNCTION 


The measure of effectiveness chosen for the 
allocation of attack aircraft was the total value of 
damage inflicted upon enemy targets. The problem is to 
maximize the total expected value of damage subject to 
certain constraints, 

(A2.1) Assume that each aircraft assigned the 1°" 
target delivers the same number of preferred weapons on 
each pass against that target. A preferred weapon 
implies that for a specific target there exists a weapon 
which is most effective in destroying that target. 

(A2.2) On any pass against a target assume that the 
target is either killed or not killed. This implies that 
a target can be destroyed on one pass. 

(A2.3) Aircraft make passes until all weapons are 
expended or until the aircraft is killed. | 

(D2.1) The base-aircraft index specifies an aircraft 
type located at a specific. base, For example, 

j = 1 denotes the set of A4 aircraft from base 


number one, 


j = 2 denotes the set of F4 aircraft from base 





number one, 


j = J - 1 denotes the set of A4 aircraft from 
base number 5, 

j = J denotes the set of F4 aircraft from base 
number S,. 
denotes the number of j°" base-aircraft 


iJ 
type assigned the (ae targeb, 1 1,2,0065 Evana 


(D2.2) N 


j = L (ee eee 


(D2.3) N, denotes the vector (Nia: Nioe coe s Nig) 
(D2 .4) N, denotes the total number of the th base- 


aircraft type available. 


(D2.5) R,, denotes the number of passes per aircraft 


th 


ar 
planned by the yon base-aircraft type assigned the i 


target. 


(D2.6) PT,» @ constant, denotes the probability 


th 


the i target is killed by exactly one pass given that 


b preferred weapons are delivered on that pass. 


.(D2.7) W, (N,) denotes the random variable which 


represents the total number of passes against the — 


target. W,(N,) will be written as W, when convenient. 


1 


i 


The range of possible values is? 


J 
W, = Wyle 25 6 ofa N 


R, 4. 
i eee 





4 be a random variable such that 
Y, =|0 if the 1°" target is not killed 


1 
1 if the 1°" target 1s killed at least once. 


(p2.8) Let Y 


P [y, = o| = EP |Y, =o|w, =w] P[w, = w] 


=-£(1-PT,, )* Pilw, =wi, 
all w ib a | 
and 


1-P[y, = 0] 


ty 
—— 
' KS 
UI 
pd 
es 
lt 


1-£(1- PT, ie 


all w b 


(D2.9) Let v, denote the pre-assigned value of the 


th 


19 target, and let V=£ v, be the total value of the 


1 


tl MH 


i=1 
target complex. 
(D2.10) Let D, denote the random variable which 
represents the value of damage to the , oh target such 


that: 
D, = |0 af Y= 8 
+ abe eae ole 
Then, 
E [D, | 3% a°P[D, = a| 
: m4P [¥, = 1] 
=v, (3 =e 1- PT, " P[W, = ‘)} 


The total value of expected damage (TD) to the target 





complex is; 


I 
(D2.11) TD = A E [ D, | 


i 
I I 
a * a ital ="? = * 
I W 
=V- LAL -- 1 - ye |W, = w | 
Z 
m= OL, 2c e's ‘ate th 


To maximize TD one must minimize the expression 


fap [ Wy = w] , which 


ries x (1 - Pr,, 


Ww 


is simply the sum of the product of the value of the Lee 


target and the probability the i°" target is not killed. 


Therefore the objective function is 
i 


(2.1) min 5. & 
i=l] 


Z(1- Pr 


all w " [" ' ai ; 


1 


1 ib 


II. THE CONSTRAINT EQUATIONS 


2 > 
pf ut A Lo Af & 
Je . f 


The eMltscactes or the objective function is subject 
to certain constraints on the resources available. The 
four constraints given below are a non-exhaustive subset 
of those which could be imposed, but are representative 
of the most important ones? 


l. The number of aircraft allocated to targets 


must be less than or equal to the total number of aircraft 








available, i.e., 


I 
(232) me k,..e W 


9 N = 0 9 
j=) 11d 


J 


where N; = the total number of yeh base-aircraft type 


= the number of j°" pase-aircraft 
th 


available, and Ny 


type assigned the i target. 


Ze The total fuel required must be less than 
or equal to the total fuel available, i.e., 
Ie 
< 
(2.3) z Pia = PF, - 


where Fy = the total fuel available for the ye base- 


aircraft type, and M4 = the fuel required for the ye 


base-aircraft type to strike the 1°" 


target. 
3. The total ordnance loading must be less than 
or equal to the ordnance available, i.e., 


I 
(2.4) £ » = B 


4=1 1 jn’ Lj Vial 


where Ban = the number of weapons of type n available for 


the jie element, and b = the number of preferred 


th 


Lapa 
weapons of type n the jen element carries to the i 


target. e 

4, There exists an upper limit on the number 
of aircraft which one is willing to lose on a given 
mission, i.e., . i 


1 





where L, = the maximum acceptable number of yeh base- 


aircraft losses, and L = the expected percentage 


i j 
attrition of the j°" base-aircraft type assigned the 
yee target. 

There are numerous ramifications to these constraints 
which will not be covered in detail here, For example 
there are additional restrictions to (2.3), (2.4), and 
to (2.5) in that, respectively, fas is constrained by 
the amount of fuel an aircraft can carry, Bate is con- 
strained by the maximum ordnance load and rack restrict- 
ions of the aircraft type, and Lys is a function of the 
aircraft type, speed, penetration altitude, the enemy 
defenses, and the number of passes made against the 
target. 


~ 


The allocation problem has now been formulated as? 


(2.1) min : Vy Tiel = PT, 4 y" p |W, = w | , 
4=1 all w 
subject to $ 

at 

(2.2) ea =N, 5 My, 70, 
I 

(293) Pas fasNay = Fy > 
ng 

(2.4) * Bs nay = Ban » 
I 

(225) & ZL Ny <L, | 
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III. NONLINEAR PROGRAMS 


Rosen ee , | 12] and Fiacco and McCormick | 20 | - 
[22 | : [ 30 | have developed algorithms which have been 
successful in solving nonlinear programming problems 
subject to linear or nonlinear constraints, 

The Sequential Unconstrained Minimization Technique 
(SUMT) of Fiacco and McCormick and the Gradient 
Projection (GP) technique of Rosen are both available 
through the IBM Share General Program Library as Share 
Distribution 3189, RAC SUMT, and Share Distribution 1399, 
SD GP 90 respectively. 

The use of SUMT is precluded in this allocation 
problem due to the non-differentiable nature of the 


objective function. That is, W, is a discrete integer- 


i 
valued random variable. The concluding remarks of the 
Share write up of RAC SUMT indicate that special sub- 
routines are being developed to handle non-differentiable 
functions. 

The program GP 90 should handle the non-differentiable 
objective function since it uses two-sided differences 
in place of the gradients. Thus it is unnecessary to 


explicitly evaluate the gradients of the objective 


function. 


11 








CHAPTER III 
THE PROBABILITY DISTRIBUTIONS 
I. ASSUMPTIONS 


The probability distributions of two of the random 
variables, the number of aircraft alive at release on 
a particular pass and the total number of passes against 
a target, are presented in this chapter. 

The following assumptions were made: 

(A3.1) An individual aircraft is assigned only one 
target per sortie, 

(A3.2) A raid is composed of one base-aircraft type, 
but more than one raid can be assigned to a target. 

(A3.3) Enemy fighters, if scrambled against a raid, 
are sent in numbers sufficient to engage each aircraft 


in that raid, 
II, DEFINITIONS 


The following definitions were useful: 

(B3..1) PD, ,(n) denotes the probability that a raid 
of size n is detected enroute from the 3°" base-aircraft 
location to the eB target. 

(UG ) PES (Non) denotes the probability that a raid 


of size n is engaged enroute from the Ga base-~aircraft 


location to the , bh target given that the raid is detected 
i 








and given that a total of N aircraft are employed in 
strike operations. 

(D3.3) PK denotes the probability that any aircraft 
in a raid is killed enroute given that the raid is 
engaged. 


(D3.4) PA, ,(n) denotes the probability that a raid 


from the ale base-aircraft location finds the 7 target 


given that n aircraft survive enroute. 


003.5) F5 (n) denotes the probability that any 


ijk 
aircraft in a raid from the j°" base-aircraft location 


survives until the aad release against the pee target 


given that n aircraft are alive commencing the first pass, 
k = 1, or given that n aircraft are alive at the (K - ier 
release for the second and subsequent passes, K = 23 


k = 1,2 5F,000 5k, ye 


(D3.6) Xs jac (Ny 4) denotes the random variable which 


represents the number of the qn base-aircraft type alive 


th 


at release on the aon pass against the i target given 


that Ny are assigned. The range of possible values is; 


Xa yy ( Ny 5) = ORV? sce For simplicity a yu! Ny 5) is 


sometimes written as X 


i)" 
i jk° 


III. DISTRIBUTION OF THE NUMBER OF AIRCRAFT 
ALIVE AT RELEASE ON THE FIRST PASS 


Under the above assumptions the event that x aircraft, 


x>0O, are alive at release on the first pass against a 


13 








target can occur in three mutually exclusive and 
exhaustive ways: 

A raid of size N is undetected enroute, finds the 
target and x aircraft among N survive target defenses 
on the first pass, 

OR a raid of size N is detected enroute, but 
unengaged, finds the target and x aircraft among N 
survive target defenses on the first pass, 

OR a raid of size N is detected enroute, engaged, 
some n aircraft survive enroute, n = x, those n find the 
target and then x among n survive target defenses on the 
first pass. 

The extension of this type of reasoning used in 
conjunction with the laws of elementary probability 
theory produced the probability distribution of the 
number of aircraft alive at release on the first pass, 
Xo42? presented in Table I. 

IV. DISTRIBUTION OF THE NUMBER OF AIRCRAFT 
ALIVE AT RELEASE ON SUBSEQUENT PASSES 

The probability distribution of Ko 4x9 k = 2, was 
obtained from the distribution of Xia by making the 
following observations: 


The stochastic process, | Xy sys k = 1s2s3peee0Ry yf 
is @ discrete parameter finite Markov chain. That is, 
a LX 5 = Xo | Xa4y = Xo Ky yo = po eee o Mpa = 3 | 


= P(X, 5, = x, | Se pel a 
14 








The number of the 3°" basge-aircraft type alive at 
release on the «eh pass against the oe target is less 
than or at most equal to the number alive on the (y-1” 
pass. That is, 


<= xX 


Xo 4x 1 jk-1 k 2 Fatty esos Ry ye 
Therefore the conditional distribution of Xin? * =r 


given that X is binomial with parameters 


le ie 
PR, y;,(m) and m as presented in Table II. 
By the theorem of total probabilities: 
N 
- —~ gid 
P [xy =z) = Bo P [x 


Bn Pp [x 


ie “S| Figce1 = n | ‘ee ™ n| 


The synthesis of these facts yielded the probability 


distribution of X 2, presented in Table III. 


13k? * > 


V. DISTRIBUTION OF THE TOTAL NUMBER 
OF PASSES AGAINST A TARGET 

The derivation of the probability distribution of 
W, necessitated one further assumptions 

(A3.5) The number of a given base-aircraft type Alive 
at release on any pass over a given target is statistically 
independent of any other baseeaircraft types alive over 
the same target on any pass. That is, Xs 4x is assumed 
independent of X,,/,,provided j # a 

Under all foregoing assumptions the probability 
distribution of Wy as presented in Table IV is derived 


in Appendix A. 
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CHAPTER IV 
THE PARAMETERS 


I, TARGET VALUE 


The value of a target may be characterized by a time 
dependent index, i.e., target value will change as the 
tactical situation changes over time, Furthermore, the 
relative values of targets in given tactical situations 
will depend strongly upon the doctrine adopted by the 
local military commander. If an index of target value 
is to provide a basis for decision, it must be formulated 
with consideration given to the whole social, political, 
economic, and ethical processes within which military 
action takes place, The intent of this section is not 
to formulate a scheme for target value indexing but to 
present some aspects of the problem which should be 
considered in any indexing plan, 

Kaplan [2] delimits areas of empirical study by 
which each element in the concept of military worth can 
be measured or appraised, These areas include values, 
objectives, welfare, achievement, and worth. 

Kaysen [3] outlines criteria for selecting target 


Systems as¢ 
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1. Military importance -- a rough classification 
of the value to enemy military operations of all types 
of equipment and supplies used by enemy forces, 

2. Percent direct (or indirect) military use -- the 
share of total output of a product or service which goes 
into military use; 

3. Depth -- an indication of the time available to 
the enemy for the organization of substitute consumption, 
alternate production, etc., before suffering military 
damage; 

4. Economic vulnerability -- includes 

a. Ratio. of capacity to output, 

b. Substitutability for processes and equipment, 
Cc. Substitutability for product, 

d. Process and plant layout vulnerability, 

e. Recuperability, 

5. Physical vulnerability; 

6. Location and size of the target system. 

In addition to the above criteria for indexing target 
value the effect on enemy morale should be considered. 

Hesse and Mitchell in ’*’Limited War Campaigns: 

Dacca Method’’ [ 10 | have derived an index of target value 
which incorporates several of the criteria mentioned 
above. They consider military targets from a balanced 
ton standpoint using weighting factors for conSumable and 


reuseable equipment, supplies, and manpower. This index 
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is quite involved, but it provides a method of assigning 
target values for most conceivable military targets. It 
Should serve well as a useful working tool for a first 


look at target value in the allocation problem. 


II, FUEL AND ORDNANCE LOADINGS 


To maximize damage to the enemy requires that air- 
craft be put over their targets with optimal ordnance 
loadings. A basic prerequisite to attainment of this 
objective is determination of fuel required to carry out 
the given mission. 

The following basic mission profile has been postulated 
to determine fuel requirements$ 

1. Warm-up and take-off, 

2. Climb to cruise altitude and cruise to descent 
point for run-in to the target, 

3. Descent to target run-in altitude, 

4, Run-in to the target, 

5. Ordnance delivery, 

6. Run-out from the target, 

7« Climb to cruise altitude and return to base, 

8. Descent to the base for landing, 

9. Landing and reserve. 

For ease of computation the following assumptions 
have been made: 


(AY.1) Aircraft carry full internal fuel and use 
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fully loaded 300 gallon external fuel tanks as required} 

(A4.2) <A fuel weight of six and one-half pounds 
per gallon was used; 

(A4.3) No fuel is burned and no distance over the 
ground is covered in descents, 

(A4.4) All-external store racks for ordnance-are 
equivalent on any particular aircraft type; 

(A4-5) There are no asymmetric load restrictions 
for launching aircraft, 

(A4-6) Five minutes fuel at normal rated power was 
allotted for warm-up and take-off, 

(AY-7) Thirty minutes fuel at sea level maximum 
endurance was allotted for landing and reserve. 

Utilizing the mission profile and the assumptions 
stated above one may formulate the basic fuel required 
equation as 

(4.1) Total fuel required = Take-off fuel + (2 x 

cruise fuel‘rate x 
(Range - Climb distance - 
Run-in distance)) + (2 x 
Run-in distance x Run-in 
fuel rate) + (Time over 
target x Target fuel rate) 
+ Landing and reserve fuel. 
After determining the total fuel required the next 


Step is to ascertain if the fuel required is less than 
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the fuel available. Successful mission performance is 
constrained by 

(4.2) Fuel available - Fuel required =0. 

Substituting for Fuel available in (4.2) one obtains 

(4.3) (Internal fuel + External fuel) - Fuel 
required =0. 

Solving for External fuel results in 

(4.4) External fuel >= Fuel required - Internal fuel. 

Once (4.4) has been solved it is possible to determine 
the number and size of external fuel tanks required and 
Subsequently to determine the number of external store 
stations available for ordnance loading. 

The maximum ordnance load may be obtained by solving 

(4.5) Ordnance load = Maximum take-off weight - 

(Basic aircraft weight + 
Internal fuel + External fuel). 

Appendix B contains a program (PROGRAM ORDLOAD) written 
in FORTRAN IV for the CDC-1604 digital computer which 
computes for a given aircraft mission the fuel required, 
the number of store stations available for ordnance, and 
the maximum ordnance loading. This program uses the 
assumptions and the algorithms outlined above for 
computations. Inputs required to the program, definitions, 
a program listing, and a sample output are given in 


Appendix B. 
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III. OTHER PARAMETERS 


The probability that an aircraft survives to make 
passes against a target is a function of the probability 
of detection by the enemy, the probability of engagement 
by the enemy, the probability of being killed enroute to 
the target by the engaging force, the probability of 
acquiring the target, and the probability of being alive 


at the weapons release point on the yeh 


pass against the 
target. Values for all of the above are required to 


obtain the distribution of W To minimize the objective 


5 a 
function one also needs values for the probability that 
the 1°" target is killed By one pass given that b 


weapons are released on that pass. These probabilities 
will be considered as input parameters, 

Derivation of each of these probability distributions 
is beyond the scope of this paper; however, the remainder 
of this section will be devoted to discussion and/or 
references to derivations of these distributions. 
Throughout the investigation of these parameters, past 
data will in some instances be adequate and reliable, and 
the parameters could be estimated by Bayesian techniques. 
World War II, Korea, and Viet Nam experiences should 
result in extensive information which may be used to 
increase confidence in the estimates used, 


Koopman /2| derives an expression for detection 
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probability both for visual and for radar detection. 
Electronic and other intelligence data coupled with 
characteristics of friendly radar should allow an 
estimate to be obtained for detection probability. Crout, 
Fay, and Harvey [5] have formulated a model describing 

a .bomber’s penetration into hostile territory in the 

face of a given area defense strategy. Specifically an 
expression, is developed to compute the probability that 
an attacking aircraft is killed before penetrating to a 
given depth. This model then accounts for the probability 
of detection, engagement, and kill. 

Simultaneously with passing through enemy area 
defenses friendly aircraft must attempt to acquire their 
assigned targets. Acquisition includes the processes 
of search, detection, identification, and flying the 
aircraft into a position to make a weapons pass against 
the target. Visual acquisition of the target is, of 
course, the most reliable and informative method of 
locating a given target. Due to circumstances such as 
weather, terrain, camouflage, aircraft speed, and aircraft 
altitude a pilot will sometimes have to rely on radar, 
infrared, microwave, or vectoring from personnel not in 
the aircraft to acquire his target. A combination of 
these methods should increase the probability an aircraft 
acquires its target. Erickson [9] presents a study on 


the visual capabilities of a pilot searching for ground 
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targets. Kamrass and Heckroth [2.5] include a section 
on target acquisition which discusses detection 
probabilities and an analysis of types of sensors and 
methods of operation which will be most effective for 
reconnaissance and strike: missions. 

Once an aircraft has acquired a target and is positioned 
for a weapons release pass it will usually be in the 
area of effectiveness of the point defenses of the target. 
These point defenses can run the gamut from no defense 
to a completely automatic system of surface-to-air 
missiles with a ground level to outer space capability. 
The probability an aircraft can successfully release 
weapons over the target is highly dependent upon the type 
of defensive weapon encountered. Selection of Aircraft 
for Tactical Air Missions [24 discusses and gives values 
for probabilities of kill and the effect of various 
penetration aids associated with several types of point 
defense weapons, 

If an aircraft survives the point defenses, it will 
release weapons against the target, The probability of 
killing the target on one pass will be a function of the 
type and number of weapons released on that pass, the 
skill of the individual pilot, the weather, the type of 
terrain, and the type of target. The Naval Ordnance Test 
Station (NOTS) has worked on the development of weapons 


and the formulation of models which allow the probabilities 
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of kill of those weapons to be computed. NOTS technical 
publications include those of a probabilistic nature -- 
Kusterer [8 | » W.B,. Simecka [ 16 | , Verry | 19 | > Kabe 
Simecka | 26 | » and Strang [ 28 | ; and those of a 
computational nature -- Weldon and Young [ 24 | » WeB. 
Simecka, et. al. | 23 | », and ’*®Conventional Air-Delivered 
Strike Weapons’? | 27 | - 

An alternative to derivation of the various 
probability distributions for all the parameters would 
be to use three estimates of each parameter, i.e., a best, 
an optimistic, and a pessimistic estimate. These will 
be referred to as BOP estimates. After obtaining 
solutions to the allocation problem using these BOP 
estimates, a sensitivity analysis should be performed 
to discover which parameters Significantly affect the 
results. Following this analysis high confidence 
estimates for the most sensitive ’’parameters’’ could 


be obtained using the distribution theory. 
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CHAPTER V 
CONCLUSIONS AND RECOMMENDATIONS 
I. CONCLUSIONS 


This analysis has demonstrated that the allocation 
of aircraft to fixed location targets is.a complex 
problem,’ If it had been reasonable to assume that 
attack aircraft make a single pass, as in strategic 
nuclear warfare, then the nonlinear model reported by 
Fiacco and McCormick [22] could have been used. However 
in limited warfare, for which this model was developed, 
it is not reasonable to assume a Single pass, Furthermore, 
we believe that a deterministic treatment of repeated 
passes would not suitably reflect the expected outcome 
in-even an approximate fashion. The substantiation of 
this assertion is dependent upon the generation of the 
probability distribution of W,° 
The objective function for this model is non- 
differentiable; this property restricts the class of 
usable and currently available nonlinear programming 


algorithms to that of Rosen [7] ° 


2g 





~ 


II. RECOMMENDATIONS 


It is recommended that: 

1. This model be coded for use on an electronic 
digital computer utilizing the algorithm of Rosen [7] ° 
2e The list of targets be grouped into type or 
character classes and an investigation made of the feas- 
ibility of using common values within each group for the 


parameters PA, ,(n); PRy yy, (1) and v, respectively. 


1 

3. The list of targets be grouped into geographical 
classes and an investigation made of the feasibility of 
using common values within each group for the parameters 
PD, ,(n) and PE, (Non) respectively. 

y, A sensitivity analysis using BOP (best, 
optimistic, pessimistic) estimates of the parameters be 
made to determime those requiring further analysis. 

5. The feasibility of invoking the Central Limit 


Theorem on the expression W, = % Z Xi 4k be determined, 


a 3 x 
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APPENDIX A 
DERIVATION OF THE PROBABILITY 
DISTRIBUTION OF W, (N, ) 


The event W, 


if and only if X 


= 0 occurs if and only if x ‘ Xian = 0, or 


14k = 0 for all j and all k. It was noted 


= xX 
i jk i jk-l’ 
for all j and all k if and only if X 


that X K = 2535hyeoesRy ys So that X — 0 


ijk 


141 7 O for all j. 


It was assumed that X, is independent of Kaya! provided 


jk 


that j # j', so that: 
P [u, = | = P(r x 0) =P[X, 4,= 0 for all j and all | 


it Xo 


= PLX, 43 


The distribution of W 


0 for all 5] =1 Mis = 0] 


‘ > 1 was obtained in a 


Slightly different manner, The event W, = 1 occurs if 
and only if for exactly one j, ERT = 1 and X42 = 0, 


but the event W, = 2 occurs if for exactly one j, 


it 
X51 = 2 and Xs 452 = 0, or if for exactly one j, Xi 41 = 1,%, y9= 
and X = 0, or if for exactly two j’s, Xo 4 = 1 and 


i j3 


X, 42 = 0. This enumeration can be continued and was up 


to Ws = 4, it produced a lengthy table from which a 


recursion relation was developed. The use of this relation 


in conjunction with the probability distribution of Xs 4x 
produced the probability distribution of W, in the following 


manner, 34 
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(D6.1) Let faco denotes for all combinations of. Where 
**combinations’’® is used in the sense of combinatorial 
analysis, e.g., Parzen [6 » then 


is |W, = 1| = z[ X51 ~ 1] 2 |X, 4p = 0 | X49 = 1 ret P(X, 3 = ) 


1 
p(w, =2j=c © Pix. - 2-2, | P(X, 42 = 8 1% y) = 2 - a, | 
j a =0 


: Pitigy, =0 | X42 ~ a, | ia P(x, 4 = 0] 


+ 2 P [X92 = 2) P [X50 = 91 X59 - i| 


Jayde 


| =1| P —_ =f i ao = 1 


PLx, na ~ 0} 


NOTE: P iis = x | Xs jKn2 = | ="O8for all x > m and 


TT 
n# jx 9 J2 


P(X, 4, 20 | ¥en - o] = 1, Table II. 
; ec P En = 3-e. a, | 


a,=0 a5=0 


p(w, = She 

J 
» «P(X, 42 =, - 22 | X41 = 3 - a) Px, 55 = a, | *je > ale a] 
» PLX, 5, =0| 4443 = a. | ah P[Xnr = 0] 


(continued on the following page) 
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By induction the general term of the probability 


distribution of W, was obtained and is presented in Table 


1 
IV, where all sums are defined in the positive sense, 


e 
i.e., the sum f£ f(y) is defined if and only if e =a 
y=d 


and the following definitions were used: 


(D6.2) RK denotes the max (R 
To 


(D6.3) x3r denotes the max (R,J) 


15) 


Summation Operators: 


rel Pax, -277 
(D6.4) B, = Tl »y 
k, =2 gt a 0 
jel a = 
(D6.5) C, = in k-ly 
k=2 ay =Jok 
Wol rel a wo ll 
(D6s6) As = rt kel, 
1 
a4=0 =? a, =0 
Wel W-a,-l aol 
(D6.7) A, = £ 3 ey ye 
ay=1 by =0 6,570 
Wel jel a oa ol a =] 
(06.8) A,= EZ  C,Jn tel, t B, J-1;, B, 
a j=j-1 t=1 b,4=0 317° 


j = 3,4 9m ieeoeand ay = W 


Dia 





Ta 
on 
k,= 


Summation Operands, 


(D6.9) G, (i) = P coer = 841 - a (1-den) - Py] 


4 [%5,2 = bey - Pp | jer ~ 7-1 - a, (i-de) - Pea | 
PIX ee x 5 =f 
ijpketl ~ Pew, 7 Penta | ijk © Menage tk, | 
[Xs 5,2 = he 1 | ioral = (Dee = Peril 
P [Xa j rt = 2 | “i447 7 Dera 
where ds; =< 0 if i # j 
1 tee 
‘ T=2 ( 
i “ae view “eel wee 
P a. = Oy - Bay | Kean = 3xK-1 7 a, | 
P La = Oa al | Ai ayr-1 7 ree 7 a4) 
P ae as | Aogir = Fry) oT P Ee 0] 
n# jy 
oe i 
(D6.11 Pp = »y TT G, (1) 
a faco t=l tm 
Jaseoosd., 
a | pul. 3= o] 
n#Jas+¢%9 by Lima 
m= 2793:5;¢43-1 ane 3 = 3: 
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(D6212) OPp= 9 eee J=2 
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Therefore under the above definitions, given values 
for R and J, the probability distribution of W, presented 
in Table. IV can be expressed in closed form for 


computational purposes, 


Ld 
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APPENDIX B 
PROGRAM ORDLOAD 


The FORTRAN IV program to compute the fuel and 
ordnance loadings is presented in this Appendix. A 
brief description of the imput data, the definition of 
variables, a program listing, and a sample output will 
be given. Data for the sample problem was selected at 
random, and the resultant output bears little semblance 
to any realistic situation that might occur. The program 
has been dimensioned to handle six aircraft type, six 
bases, and fifty targets. These dimensions may be 
expanded as necessary to suit a particular user’s need. 
Common blocks have also been established to facilitate 
incorporation of this program as a sub-routine in a 
larger allocation problem. 

The data input is presented below, Variables with 
subscript (I) require one data card for each aircraft 
type. 

Card 12 FORMAT (3110,F10.0): 

Fields 1-103 NTYPE = the number of different 
aircraft type, 1 < NTYPE « 6; 


Fields 11-202 JBASE the number of friendly 


bases, 1 = JBASE <= 65 


40 





Fields 21=-30° KTGT = the number of targets, 
1 = KTGT = 50; 
Fields 31-403 TIME = the time ( tenths of 
hours ) over the target for delivery of ordnance. 
Card 2,.0o0¢,NTYPE+1° FORMAT ( 5F10.0,110 )° 
Fields 1-10: GAXTOW(I) = maximum gross aircraft 
take-off weight (pounds) $ 
Fields 11-203 BACFTW(I) = basic aircraft weight 
(pounds) $ 


Fields 21-30% ENTFW(I) maximum internal fuel 


ti 


Capacity (pounds) > 
Fields 31-40 EXTFW(I) = maximum external fuel 
capacity (pounds) § 
Fields 41-503 GAXBMW(I) = maximum gross ordnance 
load (pounds) $ 
Fields 51-603 IRACK(I) = maximum number of 
external store racks on an aircraft; 
Card NTYPE+2,..0,(2xXNTYPE)+1°2 FORMAT (5F10.0)° 
Fields 1-103 FTAKE(I) = fuel (pounds) required 
for warm-up and take-off; 
~s Fields 11-203 FCLIMB(I) = fuel (pounds) 
required for climb to cruise altitude; 
Fields 21-303 FLOIN(I) = loiter fuel rate 
(pounds per minute)== not required in the current 


formulation of the program; 
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Fields 31-40% DASHSP(I) = speed (knots) over 
the target at military power; 
Fields 41-50: FLAND(I) = landing and reserve 
fuel (pounds). 
Card (2xNTYPE)+2,...,(3xNTYPE)+1:2 FORMAT (2F10.0): 
Fields 1-103 DISTML(I) = distance (nautical 
miles) covered in the climb to cruising altitudes 
Fields 11-20: DISTRI(I) = run-in distance 
(nautical miles) to the target. 
Card (3NTYPE)+2: FORMAT ( 2F10.0): 
Fields 1-10: ALT1 = cruise altitude (thousands 
of feet); 
Fields 11-203 ALT2 = run-in altitude (thousands 
of feet); 
Card (3xNTYPE)+3,...,(4xNTYPE)+223 FORMAT (4F10.0)? 
Fields 1-102 SFC(I).= sea level cruise 
specific fuel consumption (pounds per nautical mile = 
ppnm ) 5 id 
Fields 11-203 SFCMIL(I) = sea hevel military 
power specific fuel consumption (ppnm); 
Fields 21-302 SFCK(I) = cruise specific fuel 
consumption altitude correction factor (ppnm per thousand 


feet altitude); 
Fields 31-40: SFCMILK(I) = military power 


specific fuel consumption altitude correction factor 


(ppnm per thousand feet altitude). 
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Card (4xNTYPE)+3,..0,5 (4xNTYPE3+2+(JBASEKKTGT)$: FORMAT (F10.0); 
Field 1-103 RANGE(J,K) = range (nautical miles) 
from Jth base to the Kth target, 


Definitions of other variables used in the program 


ares 

FCRUSA(I) = cruise altitude fuel consumption rate 
(ppnm) ; 

FCRUSL(I) = run-in altitude fuel consumption rate 
(ppnm) ; 


FDASHM(I) = run-in altitude military power fuel 
consumption rate (ppnm)>5 

FAVAIL(I) = maximum total fuel available for 
aircraft type I, 

FRQRD(I,J,K) = fuel required for aircraft type I 
to complete mission from base J to target K and return, 

EXTFRQ(I,J,K) = external fuel required for aire 
craft type to complete mission from base J to target K 
and return; 

KRACK(I,J,K) = the number of external store racks 
available on aircraft type I to target K from base Jj 

ORDLD(I,J,K) = ordnance load (pounds) aircraft 
I can carry to target K from base J. 

A complete program listing and a sample output, which 

uses the data at the end of the program listing, are given 


on the last four pages of this Appendix. The sample 
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output tables the aircraft number, the base number, 

the target number, the range from base J to target K, 

the total fuel required (lbs.), the internal dnd txteowa 
fuel required (lbs.), the number of ordnance racks available, 


and the ordnance load capability (pounds). 
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